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Abstract 

We propose a simple laboratory experiment to measure the speed 
^ ■ of propagation of gravitational phenomena. We consider two masses 

placed at two different distances from a detector of gravitational force. 
They are placed so that the static gravitational force they produce at 
^ ■ the detector exactly cancels. We show that if the masses are made 

. to oscillate, then the force on the detector no longer vanishes and is 

proportional to the relative time delay it takes for the propagation of 
I the gravitational changes from the masses to the detector. This time 

■ delay is inversely proportional to the speed of gravity. A measurement 
, of the finiteness of this speed would be the first confirmation of one 

dynamical aspect of Einstein's theory of general relativity, that gravi- 
tational effects are not instantaneous and propagate at a finite velocity 
and a direct confirmation of the fundamental notion that there is no 
^ I action at a distance. It is fully expected that the speed of gravity is 

■ equal to the speed of light. 
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1 Introduction 

Static gravitational phenomena has been described by the theory of New- 
ton since the 17th century with striking accord with observations in the 
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solar system and in the laboratory. In the early 20th century this theory 
was superseded by the dynamical theory of Einstein, which is usually called 
general relativity. General relativity is considered one of the summits of 
human intellectual advances and revolutionizes our understanding of the na- 
ture of gravitation and its interaction with matter, and of the nature of space 
and time. General relativity was successful in explaining certain kinemati- 
cal discrepancies between Newton's theory and observations, especially the 
bending of starlight by gravitating bodies and the surfeit in the advance 
of the perihelion of Mercury. Up to now, direct verification of the theory 
of general relativity has been restricted to its kinematical predictions, that 
matter causes the curvature of space-time and material bodies move along 
geodesies of the curved space-time manifold. However, general relativity also 
predicts the complete dynamical aspects of gravitation. Foremost of these is 
the prediction that gravitation propagates through space-time in the form of 
quadrupolar gravitational waves, which travel at the same speed as the speed 
of light. There has been indirect evidence for the existence of gravitational 
waves through the observation of the decay of the orbit of co-orbiting binary 
pulsars prj. The rate of loss of energy in the system matches exactly with the 
prediction of how much energy is taken away through quadrupolar gravita- 
tional radiation. However, these waves have not yet been directly detected. 
There is a prodigious effort underway to construct enormous, very expensive 
gravitational wave detectors, which would be able to directly detect gravita- 
tional waves through observing their effects on the contraction and expansion 
of the lengths of the arms of massive interferometers |2j. 

In this article, we propose an experiment, which as a consequence of the 
expected finite speed of propagation of changes in the gravitational field, 
allows us to measure this speed. We suggest that the speed of gravitational 
phenomena can be measured through table top experiments. Our analysis 
depends on the idea that if one is able to detect the static gravitational 
field of a massive body in the laboratory, then moving that body induces 
detectable changes in the gravitational field that propagate to the detector 
with a finite delay. 

Suppose we have a detector, which reacts to the gravitational field. If we 
change the gravitational field by moving masses in the neighbourhood of the 
detector, the detector will register the change. However, the change detected 
at the detector occurs after a certain delay, that equals to the travel time for 
the gravitational change to propagate from the mass to the detector. This 
delay is equal to the distance between the body and the detector divided by 
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the speed of propagation. Moving just one body in the presence of a detec- 
tor simply introduces the corresponding changes in the gravitational field at 
the position of the detector at the appropriately advanced time. However, 
with two or more bodies moving in the presence of the detector, it is easy 
to conceive of a situation where the changes in the gravitational field at the 
position of the detector will contribute constructively or destructively. The 
interference will depend on the relative distance of the two bodies from the 
detector. Observing this interference should be possible and would allow for 
the determination of the speed of propagation of the changes in the gravi- 
tational field based on the model of general relativity, as we will explain in 
the rest of this article. Throughout this article we will use the symbol Vg 
to denote the speed of gravity, expecting full well that Vg — c, the speed of 
light. 



2 Theoretical analysis 

Consider the following geometry. A stationary detector of mass m at position 
X is separated from a moving point source of mass M at f{t) , hence they are 
separated by a distance r = \x — f{t)\. The gravitational force on the detector 
is given by Newton's law as 

|F| = . (1) 

However this expression is not actually correct. In fact, the force at position 
X and at time t corresponds to the Newton expression but with the source 
mass at the retarded position f(t'), where t' — t — \x — f{t')\/vg. We will 
give a proof of this in the next section, but first we will extract some useful 
consequences assuming this delay. 

If just one mass moves in the presence of the detector, the effects are felt 
at the detector with an overall delay, and it is difficult to conceive of how to 
measure this delay. If we arrange two masses, Mj, placed along the x axis, 
at distances rj(t), for i = 1, 2, on either sides of the detector, which is placed 
at the origin, then the x component of the force on the detector is given by 

F(xt)-^^-^^ (2) 

where t[ = t — riit'^jvg. Taking rj(t) = i?j -|- 5^ sin cut = + sin cut) and 
assuming \ei\ -C 1, then = t — Ri{\^ei%m.ujt'^lvg ~ t — Ri/vg-\-o{ei). This 
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yields 



F,{x,t) 



+ 



+ ei sinu;(t;)))2 ~ (i?2(l + ea sina;(t - i^aM)))' 

GmM2 GmMi 
+ — 52 — 2e2sma;(t - i?2/fg) 52 — 2ei sma;(t - (3) 

-^2 -Ki 

to first order in e^. Let us first analyze the most evident case, '^'^^^^ = 
= ^ and ei = e2 = e. Then we find 

GmM , , , , , , 

F^[x,t) = 2e(sincj(t - i^sM) - smuj{t - i^i/t;^)) (4) 

to first order in e. Using the simple trigonometric identity, sin a — sin b = 
2 sin ^ cos gives 

i^x(x,t) = ^^4esin(a;(/?i-i?2)M)cos(u;(t-(/?i + i?2)/t^5))- (5) 

This expression is worth emphasizing. The relative delay. At = — 
R2) I Vg tempers the forced oscillations induced by the cosine term. There is an 
overall delay, {Ri + R2)/vg, but it is difficult to see or squarely unobservable. 
Putting some physical numbers into the formula, for masses of the order of 
kilograms and distances of the order of meters, the factor in front, ^^^^ ~ 
lO^^^Newtons. Ahhough very small, it was perfectly measurable even in the 
time of Cavendish. We now have a further attenuation of the amplitude of the 
force by e sm{u{Ri — R2)/vg). The force is periodically varying with temporal 
dependence cos(a;(t — + R2)/vg)). If we imagine e ~ 1/10, a; ~ 10^ and 
At ^ 10"^ we find esin(w(i?i - R2)/vg) ^ 10"^ Thus, the effect will be 
observable if the detector corresponds to an oscillator of resonant frequency 
w, and of sufficiently high Q value. Resonant mechanical detectors, such as 
torsion balances and micro-cantilevers are known to have Q values of the 
order of 10^ [31, ?]. Therefore, with specifically designed experiments, it does 
not seem to be out of the question that the effect could be observable. 

Taking the general case, = F± 5F/2 and = e ± 5e/2, with the ± 

signs correlated with i = 1 or 2, and allowing for a possible small, intrinsic 
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asynchronicity in the motion of the two masses 5t, we get (advancing the 
time argument and expanding to first order in uj{At — 5t)) 

F,{x,t+{Ri + R2)/vg) = 
= SF + {Fe + SF5e/4:)2smuj{{Ri - R2)/vg - 5t) cos{cot) 

-{5Fe + F6e) cosa;((i?i - R2)/vg - 6t) sm{ut) 
^ 6F + {Fe + 6FSe/A)2uj{At - 6t) cos(wt) - {6Fe + F6e) sin(wt) 
= 6F + {A^ + B^y/^cos{ujt + a). (6) 

Here A = {Fe + SFSe/4:)2uj{At-St) and B = {5Fe + F5e) and tana = B/A. 

This expression is a function of three parameters, 6F, 6e and 6t. These 
parameters should be thought of as adjustable experimental parameters that 
can be brought to zero by making changes in the experimental apparatus. 
To make 6F = and 6e = 0, we simply have to adjust two parameters, 
the equilibrium position and the amplitude of the oscillations of one of the 
masses, until the amplitude of the driving force, and consequently that of 
the resonating detector, becomes minimal. Experimentally, this requires a 
search in a two dimensional parameter space. When attained, we have B = 
0, a = and A = Fe2co{At — 5t). It is further possible to adjust the 
intrinsic asynchroncity so that 5t = At which makes the entire driving force 
vanish. But this procedure then just gives us a measure of At. We can 
then directly extract the speed of gravity by measuring Ri and i?2- Indeed, 



Assuming the expected relativistic invariance of the propagation of gravita- 
tional effects, the field at the detector can be calculated from the gravitational 
analog of the Lienard-Wiechert [6] potentials for point sources. Indeed, the 
empty space, matter free solution to the Einstein equations for the gravita- 
tional metric field is given by the Minkowski metric 



vg = - R2)/At. 



3 Mathematical justification 



/ 1 \ 
0-10 

0-10 

\ -1 / 



(7) 



and adding a small, perturbing source yields 



(8) 
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We get in first approximation, in the harmonic gauge [7] 



UK. = ^5.. (9) 

and Sf,^ = T^^ - {l/2)T\ri^^, and □ = d'^/dt'^ - v^gV ■ V with all derivatives 
taken with respect to the fiat coordinates. Vg is the speed of gravity, which 
is of course expected to be the same as the speed of light. 

The solution to this equation is obtained by using the retarded Green's 
function 

Ku{x,t) = ^ d^x'dt'S^,{x\t')^ ,^ ' " (10) 

Vg J \X—X'\ 

where the delta function imposes the retarded time condition, viz. the change 
in the metric, /i^i,, at position x at time t is given by Newton's law for the 
point mass at position x' and at time t' where t' = t — \x — x'\/vg. Clearly 
\x — x'\/vg represents the time it takes for the signal to travel from x' to x. 

However, we will be considering a mass of finite size, moving rigidly. 
Suppose a point in the mass moves according to the trajectory r(t), then all 
other points in the body are displaced with respect to this point by a time 
independent displacement, r. The energy momentum tensor for a point mass 
fh is given by 

TA^.t) = ^^^^§^S%S-m) ^ m(3,{t)^,{t)6-\x-r{t)) (H) 

where the last approximation is valid in the non relativistic limit, = 
{l,Vi/vg). Then 3,,^ = m{(5^{t)[5y{t) - {l/2)r]f,^)S^{x - f{t)), replacing m 
d^rp, r{t) — )■ r{t) + r and integrating over d^r, we get 

h,u{x,t) = j d^rp{f) j dt'{f3,{t')f3,{t')-il/2)r],,)x (12) 

^ f - ^(^') - ^^(^ -t'-\x- X'\/Vg) ^^g^ 

/ \x — x'\ 



-AG 



d'rp{r) / dt'{P^{t')P,{t') - (l/2)r/^,) x (14) 



6{t-t'-\x-r{t')-f\/vg) 
^ |f _ f(f ) - f\ • ^^^^ 
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Now the assumption that we must make is clear, if |r| <^ |x — r{t')\ then we 
can neglect r to lowest approximation. Writing M = J d^rp{f), we find 

Vg J \x r[t ) I 



(16) 

Then for example, we find, dropping all velocity dependent terms, 

-2GM 

hoo{x, t) = — — — — — - = hii{x, t) = h2{x, t) = /i33(f , t) (17) 
Vg\x — r[r)\ 

where t' is the retarded time satisfying t' = t — \x — f{t')\, which is the lowest 
order approximation to the Schwarzschild-like metric 

This means that the correct potential is simply the usual Newtonian potential 
but corresponding to the source at the retarded position, as we have used 
above. 



4 Perspectives for measurements 

Here we analyze the perspectives for the possibilities of observing the effect 
and determining the speed of gravity. 

4.1 Torsion balances, cantilevers and other mechanical 
detectors 

Torsion balances are the oldest, most venerable detectors of gravitational 
phenomena [3]. The greatest problem, however, is that the amplitude of 
the swing of the torsion balance can be appreciable, in fact non-negligible 
compared to the amplitude of the oscillating masses. In this case, one must 
take into account the movement of the detector, which we have assumed to 
be negligible. Micro-cantilevers |1],[5] on the other hand, have very little 
intrinsic movement. Superconducting gravimeters |8] may well be better 
suited. Such devices levitate a superconducting test mass in the magnetic 
field of persistence currents of superconducting loops. As the currents are 
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exactly stable, because superconducting persistence currents do not decay, 
the levitating force has no variation, and hence the system is ideally suited 
to detecting changes in the gravitational field. Such devices could be ideally 
suited for detecting our proposed effect. 

In the end, the Q value of the resonant detector will govern whether or 
not it will be able to detect the effect. Q values of mechanical detectors lie 
in the range 10^ to 10^ which is at the limit of what is required. We have 
a driving force with an amplitude that is approximately 10~^ smaller than 
what is an easily detectable gravitational force. However, we imagine our 
masses made of the most dense materials available, a kilogram of Osmium 
corresponds to a sphere of radius 3.25 centimetres. Gold is not that far off. 
Then the distance to the detector could be cut to about 30 centimetres while 
the oscillation amplitude could be taken to be about 3cm. Then the force 
is increased by a factor of 10. Thus, specialized detectors and experimental 
set up seem to indicate that detecting our effect is not out of the realm of 
possibility. 

4.2 Gravitational wave interferometers 

The very large arm, interferometric gravitational wave detectors [2] could 
be conceivably used to detect our predicted time delay. Indeed, at first 
sight, they seem ideally suited. The experimental idea would be to bring 
two masses in the neighbourhood of one of the arms and oscillate them 
accordingly to obtain a time dependent delay in the light travel time in that 
arm due to the Shapiro time delay p5]. However, there are two problems to 
this implementation, one of which is surmountable, while the other does not 
seem to be. 

4.2.1 Null condition 

The first problem has to do with the fact that the effect that we wish to 
detect is encoded in the metric. The Shapiro time delay occurs because 
the null geodesies for light, in the presence of a mass are lengthened due to 
the bending of light. But the null geodesies are properties of the metric of 
the theory. The metric of the theory contains the Newtonian gravitational 
potential, and not the gravitational force. Bringing two masses on either side 
of the trajectory of a light ray, the gravitational potentials at the position of 
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the light ray simply add. The gravitational potential is given by 



GM 



(19) 



r 



which is negative definite. Therefore, the null condition is not easily obtained 
for the potential. In contradistinction, if the detector is sensitive to force the 
null condition is easy to establish, the forces cancel being in opposite direc- 
tions. It would also be possible for the potential, if negative mass existed, 
the potential of a negative mass has the opposite sign from that of a positive 
mass. 

However, it is easy to conceive of an effective negative mass. The potential 
of a spherical body of uniform density, with a smaller spherical cavity in 
it, is exactly as if the volume of the cavity was simultaneously filled with 
equal, positive and negative mass density. Thus the mass density effectively 
corresponds to a full, solid sphere of positive mass density, and a smaller 
sphere located at the position of the cavity filled with negative mass density. 
If we write Mq as the positive mass of the sphere without the cavity, and 
—Ml the mass of the cavity as if it were filled with negative mass, and tq the 
position of the centre of the cavity, the potential at a position x is given by 



This is simply the linear superposition of the potential of the full, cavity 
free, larger spherical body of positive mass, plus the potential of the smaller 
spherical cavity but now filled with negative mass. 

Now imagine rotating the spherical body on an axis that passes through 
its centre, which is a distance Rq from the detector, but an axis that is 
perpendicular to the separation direction to the detector. Such a source 
would produce a time dependent negative mass potential at the position of 
the detector. The total potential at the position of the detector, adding in 
the potential of an ordinary sphere of mass M2 which is moving in a mirror 
reflection of the motion of the cavity, using the conventions from Equation 




(20) 
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([3]), is given by 



GMo GMi 
$ = + 



^0 i?i ((1 + ei sin(wt;))2 + el cos^{ut[))f^^^ 
GM2 

i?2((l + 62 sin(u;t'2))2 + el cos2(a;t^)))(V2) 
+ 



i?o i?i(l + 2eisin(u;t;) +e2)(V2) 



i?2(l + 2e2 sin(u;t'2) + e2)(i/2)- 



(21) 



Interestingly, the fluctuations in the transverse direction only contribute at 
second order, and in fact in a time independent manner. Expanding to first 
order we find 

$ ^ -^ + -^(i + e,smiut[) + o{el)) 



i?2 



l + e2sin(u;4) + o(e^)) (22) 



Now taking Mi/Ri = M2/R2 = M/R, ei = €2 = e, creates exactly the null 
condition for the time dependent part of the gravitational potential at the 
position of the detector, and thus solving this proble nfl 

GMn GMe , , , s //XX / s 

$ ^ -— ^ + ^(sin(a;t;) - sin(a;t'2))- (23) 

Also, it should be pointed out that this expression is strictly the Newtonian 
potential. We do realize that rotating extended bodies does not correspond 
to rigid body translations. For example, the fact that the larger spherical 
body is rotating does affect the energy momentum tensor of the matter fields 
and the consequent metric. Additionally, rotating the cavity and the mirror 
reflected mass do not correspond to rigid translations of these masses. The 
full metric tensor in the Einsteinian theory contains additional terms which 
are due to the rotating, though stationary, energy-momentum distributions. 



^We thank G. Kunstatter for pointing out that it is also possible to obtain the null 
condition by taking ordinary masses but simply changing the phase relation of the oscil- 
lation by TT. Specifically if Mi was just an ordinary mass (not a cavity) this changes the 
sign of the second term on the RHS of Equation (|22|) . Then we must redefine Mq and i?o, 
and if we take ei = —£2 = e, we achieve the null condition. 
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The metric will have Kerr (TU] geometric type corrections, and the dynamics 
as given by the Newtonian potential, are only an approximation. These 
corrections are, in any case, very small as the velocities involved are non- 
relativistic. 



4.2.2 Metric, the dimensionless potential and the calculation of 
the Shapiro time delay 

However, the second problem is that the effect is unfortunately undetectably 
small, and this problem seems insurmountable. It is not the simple New- 
tonian gravitational potential $ that appears in the metric, but ^/Vg. The 
potential for kilogram sized masses and meter sized distance is of the order of 
10~^^ Joules/kilogram = 10~^^m^/s^, which is perhaps measurable, but what 
appears in the metric is the dimensionless potential obtained by dividing by 



V 
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lO^^m^/s^ which implies the potential in the metric is in fact, of the 
order of ~ 10"^®. Thus, the Shapiro time delay corresponds to wavelength 
dephasing by MG/Vg ~ 10~^®m which is not measurable. Specifically, the 
metric at x,t for a mass M at position f{t') [t' the retarded time), in our 
approximation and in cartesian coordinates, is given by 

(24) 

There will be an additional contribution from the spherical cavity with ef- 
fectively negative mass, which will simply subtract from that of the positive 
mass, and hence we shall not treat explicitly. 

Consider the arm of the interferometer along the and the separation 

of the mass from the interferometer along the y axis. Nothing depends on 
z, and it can be dropped. Thus x = xx and r{t') = r{t')y and \x — r{t')\ = 
(x^ _l_ r(t'YY^'^. The light travels on null geodesies, dr^ = 0, which gives 

^GM \ 2 . 2 _ A 2GM \ 2 

- ^2(2,2 +^(^,)2)l/2 ) - ^1 + ^;2(a;2+^(^/)2)l/2 J ^25) 

that is (in lowest approximation) 

/ GM \ , 
Vndt = 1 H — -— ■— , dx. (26) 
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The first term gives the usual Euclidean result while the second term gives 
the additional time that it takes - the Shapiro time delay - for the light to 
travel over the distance dx at the space-time position x, t. We should not 
forget that t = t{x), (or vice versa) implicitly a function of x defined by 
Equation fl26|) . 

In the experimental setting, if the oscillating masses are situated at Li 
and L2 from the semi-refiecting mirror and the end mirror, respectively, of 
one arm of the interferometer, we must integrate over x from to Li, add to 
it the integration from to L2 and then multiply by 2, to get the total light 
travel time. We assume r(t') = ro(l + esincut'), and we take e — )■ 0. Then we 
get t' ^ t — (x^ + r^y/"^ /vg + o(e), and we have the time delay, on top of the 
Euclidean value VgT^ = L, (in obvious notation and to be taken first order 
in e) 

^ GM 

Jo t;2(a;2 + r2(l + esinu;(t-(x2 + r2)VVt;,))2)i/2 ^ 

^GAf f 1 (rlesinujit-ix' + riy/yvg))',^^^ 
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The t appearing in Equation (1271) is a function of x, but since the delay is 
already first order, we can use the zero order expression t{x) = x/vg. There is 
no point in doing this further calculation. To get an estimate of the order of 
the size of the RHS in Equation (j27l) we replace the sin ^ 1, assume tq <^ L 
and use simple dimensional analysis replacing the integration by a simple 
factor L to get 

„ ^'^(l-t) (28) 



L2 



The length L is almost 11km for the case of Ligo [2], thus the second term(and 
in actual fact the time dependent part) is incompletely negligible. Thus we 
find VgATi ^ MG/Vg ~ 10~^^m which is simply too small to be measurable. 

Nevertheless, it may still be practical to use these interferometers by 
making the mirrors oscillate due to the moving external masses. Indeed, the 
end mirror is so highly isolated from any mechanical vibrations, that even 
quantum mechanical uncertainty in its position has to be taken into account, 
|llj . Thus oscillations induced in the end mirror directly by gravitational 
force might be easily measurable. 
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4.3 Quantum resonance detection, quantum interfero- 
metric devices 

Such devices are sensitive to the gravitational potential We mention that 
they have exceptionally large Q values, of the order of 10^° [12]. We will men- 
tion only one such system, that corresponding to the experiment q-Bounce 
[T3], [H]. Here ultra-cold neutrons are trapped in a container. They are 
cold enough so that they do not penetrate the walls of the container. They 
only feel the gravitational potential of the earth. The quantum mechanical 
wave-functions in a linear gravitational potential are the Airy functions, and 
implementing the boundary conditions gives discrete energy levels that are 
separated by pico-eV. The experiment has already succeeded in observing 
resonant transitions between the energy levels. Pico-eV corresponds to fre- 
quencies of a few hundreds of Hertz. This is ideally suited for physically 
oscillating macroscopic masses. Indeed, if the potential was made to os- 
cillate at the resonant frequency of the transition, it could be possible to 
induce transitions between the neutron energy levels. The gravitational po- 
tential of the earth is approximately 10 m/sec^. The gravitational potential 
of kilogram sized masses at distances of the order of meters is approximately 
lO^^^m/sec^. Thus the Q value of the system must be in the same range, 
10^^ ~ 10^^ in order for the effect to be observable. 

5 Conclusion 

We have proposed a simple experiment, which would in principle allow the 
determination of the speed of propagation of gravity. It would be an impor- 
tant, confirmation and validation of Einstein's general relativity if one could 
observe that gravitational disturbances are not instantaneously transmitted 
throughout all of space but propagate through it at a finite velocity, and 
indeed, to measure that velocity. 
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